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Due to Earth’s revolution around the Sun, the expected scattering rate in direct dark matter 
searches is annually modulated. This modulation is expected to differ between experiments when 
given as a function of recoil energy Er, e.g. due to the gravitational focusing effect of the Sun. A 
better variable to compare results among experiments employing different targets is the minimum 
speed Vmin a dark matter particle must have to impart a recoil energy Er to a target nucleus. It is 
widely believed that the modulation expressed as a function of Wmin is common to all experiments, 
irrespective of the dark matter distribution. We point out that the annual modulation as a function of 
Vmin, and in particular the times at which the rate is maximum and minimum, could be very different 
depending on the detector material. This would be an indication of a scattering cross section with 
non-factorizable velocity and target material dependence. Observing an annual modulation with at 
least two different target elements would be necessary to identify this type of cross section. 


INTRODUCTION 

Dark matter (DM) is the most abundant form of mat¬ 
ter in the Universe and its nature still remains a mys¬ 
tery. More than 80% of the mass of our galaxy resides 
in a spheroidal DM halo, which extends well beyond the 
visible disk. Efforts to detect new elementary particles 
which could constitute the DM are multi-pronged. 

Direct DM detection experiments attempt to detect 
the energy deposited by DM particles in the dark halo 
of our galaxy when they collide with nuclei inside a de¬ 
tector. An unmistakable signature of the expected DM 
signal is an annual modulation of the rate caused by the 
rotation of Earth around the Sun [1]. For DM velocity 
distributions that are smooth and isotropic in the galactic 
frame at Earth’s location, the expected differential rate 
for DM scattering onto a target nuclide T in all direct 
DM detection experiments could be well represented by 
the first two terms of a harmonic expansion (see e.g. [2]), 


rate among experiments employing different targets. A 
better variable is Umin, the minimum speed a DM particle 
must have in Earth’s rest frame to impart a recoil energy 
Er onto a target nucleus. It is typically assumed that 
tmax and tmin as functions of v^in do not depend on the 
target, and consequently they can be used to test the 
agreement between putative DM signals across multiple 
detectors. 

Here we point out that, in general, the annual modula¬ 
tion of the rate as a function of Umin can vary significantly 
for different target materials. Specifically, we show that if 
the velocity and target dependence cannot be factored in 
the differential scattering cross section, observables asso¬ 
ciated with the modulation, such as tmax and tmin, niay 
be highly target dependent. Our observation does not 
rely on any assumption regarding the DM distribution. 
As an illustration, we show that for DM particles with 
a magnetic dipole moment tmax and tmin depend on the 
target material. 


J^(ER,t) = SoiEn) + 5m(ER) cosf-^(t - to)) . 
dE-R V1 year J 

(1) 

Here Ur is the nuclear recoil energy and to is the time 
at which the speed of Earth with respect to the galaxy 
is maximum, close to June 1®*. At high Er, with Am 
positive to equals the time tmax at which the rate is max¬ 
imum, while tmin, the time at which the rate is minimum, 
is six months apart from tmax (except for a shift of about 
a day due to the eccentricity of Earth’s orbit). At low 
Ar, Am could become negative, implying to equals tmin 
instead of tmax (see e.g. Fig. 8.2 of [3]). Anisotropies in 
the local DM velocity distribution modify this picture, 
in particular by making tmax and tmin energy dependent. 
The gravitational focusing (GF) of DM particles due to 
the Sun inherently makes the local DM halo anisotropic 
[4]. Ref. [5] has shown GF to have a significant effect on 
the phase of the modulation at low enough recoil energy. 

Since Ar depends on the target nuclide mass, it is not 
a good variable to compare the annual modulation of the 


DM SIGNAL AND ITS MODULATION 


For the spin-independent and spin-dependent contact 
interactions usually considered, the differential scattering 
cross section is 


dap 

dE^ 
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with rriT the target nuclide mass, fj-T the DM-nucleus 
reduced mass, gt the total cross section for a point-like 
nucleus, and Ft{Er) the appropriate nuclear form factor. 
The differential scattering rate per unit target mass. 
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with p and m the local DM particle density and mass, 
respectively, and Ct the nuclide mass fraction in the de¬ 
tector. Here we defined the velocity integral 



where f{v,t) is the DM velocity distribution in Earth’s 
frame. The time dependence arises due to Earth’s rev¬ 
olution around the Sun. The modulation of the rate 
in Eq. (4) is determined by the time dependence of 
r](vuiin,t), which is common to all experiments. There¬ 
fore, for the interaction in Eq. (2), tmax and t„iin for 
fixed Umin do not depend on the target material. This 
remains true for other differential cross sections where 
the velocity and target dependences can be factored. In 
general, however, the differential cross section can con¬ 
sist of multiple terms with different velocity dependences 
and target-dependent coefficients, e.g. with DM parti¬ 


cles interacting through a magnetic dipole [6-31] or an 
anapole moment [6, 28-35]. It also happens with some 
of the interactions described by the effective operators 
studied e.g. in [36-41] (see [26, 42-45] for explicit for¬ 
mulas of scattering amplitudes). In this case the annual 
modulation of the rate can be strongly target element 
dependent. 


AN EXAMPLE: MAGNETIC DIPOLE DM 

Here we study in detail the case of a Dirac fermion 
DM candidate x that interacts with nuclei through a 
magnetic dipole moment with interaction Lagrangian 
■ The differential cross section for 
elastic scattering off a target nucleus T with Zt protons 
and spin St is 
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with a = e^/47r the electromagnetic constant, rrip the 
proton mass, At the nuclear magnetic moment in units 
of the nuclear magneton ej{2mp) = 0.16 GeV“^, and 
^'R(^^min) = ■ The first term is due to DM 

dipole-nuclear charge interaction, and the correspond¬ 
ing charge form factor coincides with the usual spin- 
independent nuclear form factor while the 

second term is due to the dipole-dipole interaction and 
has a nuclear magnetic form factor Em,t(^r) (both form 
factors are normalized to 1 at zero momentum transfer). 
We compute the cross section with the formalism and the 
form factors provided in [42, 43]. 

The differential cross section in Eq. (6) contains two 
terms with different velocity dependence: one with the 
usual \/v^ factor and another independent of v. The 
differential rate (see Eq. (3)) is thus also a sum of two 
terms, one containing r]{vniin,t) in Eq. (5) and the other 
containing 


b(nmin5 i') — 



V f{v, t) d^T . 


(7) 


For purposes of illustration we assume the Standard Halo 
Model (SHM), in which the DM velocity distribution is 
an isotropic Maxwellian on average at rest with respect 
to the galaxy (see e.g. [46] for details). Under this as¬ 
sumption the two velocity integrals p and fj have a very 
different time dependence. This can be seen in Fig. 1 
where their time of maximum Tmax and minimum Tmin 
are shown. Instead of Tmin, we plot — Tmin where 


Lnin is the time six months apart from Tmax- Fig. 1 
shows the effect of including (solid lines) and neglecting 
(dashed lines) GF. Neglecting GF, Tmin is almost indis¬ 
tinguishable from f„iin, and thus is not shown. Unless 
otherwise stated, we include GF and the eccentricity of 
Earth’s orbit in our calculations. Notice that Tmax (Anin) 
as a function of Umin coincides with the maximum (mini¬ 
mum) of the differential rate, tmax (^min)) only when the 
velocity and target dependence can be factored in the 
differential scattering cross section. 

The modulation of the differential rate depends on the 
interplay of the terms containing p and p. Since the 
relative coefficients are in general target dependent, as 
well as DM particle mass dependent, the modulation also 
depends on the target and on m. Let us denote with r and 
f the terms of the expected differential rate containing p 
and p, so that dRT/dER = r + r. Fig. 2 shows the rate 
fractions / = r/(r + r) and f = f/{r + f) as functions of 
^^min for four different target elements (fluorine, iodine, 
xenon, and germanium) employed by current DM direct 
detection experiments. For target elements with more 
than one isotope (Xe, Ge), we sum Eq. (3) over isotopic 
composition. Solid (dashed) lines in Fig. 2 correspond 
to a 100 GeV (1 TeV) DM particle. Notice that because 
of the negative sign in one of the dipole-charge terms 
in Eq. (6), r and / are allowed to take negative values. 
When this happens, / > 1 since / -I- / = 1. 

Figs. 1 and 2 can be used in combination to understand 
the target-dependent behavior of the time of maximum 










3 



0 100 200 300 400 

''min [km/s] 


FIG. 1. Time of maximum Tmax (top) and minimum Tmin 
(bottom) of rj and 77 in the SHM, as functions of Umin, including 
(solid lines) and neglecting (dashed lines) GF. The bottom 
panel shows Tmin — Tmin, with fmin the time six month apart 
from Tmax- Neglecting GF, Tmin is almost indistinguishable 
from Tmin, and thus is not shown. 
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FIG. 2. Rate fractions / = r/(r + f) and / = f/(r + f) for 
fluorine, iodine, xenon, and germanium. Solid (dashed) lines 
for m = 100 GeV (1 TeV). 


imax and minimum of the rate for magnetic DM, 
shown in Fig. 3 for scattering off fluorine, sodium, iodine, 
xenon and germanium. Solid (dashed) lines correspond 
to TO = 100 GeV (1 TeV). Also shown In Fig. 3 are the 
Ar thresholds for LUX [47] (3.1 keV, employing Xe), 
SuperCDMS [48] (1.6 keV, Ge), DAMA [49] (6.7 keV 
for Na and 22.2 keV for I), and PICO [50] (3.2 keV, 
F), translated into Umin for tot averaged over isotopic 


composition and elastic scattering for to = 100 GeV. For 
larger to, these thresholds move to lower Umin values. 

Fig. 3 shows that tmax and tmin become essentially tar¬ 
get independent above Umin — 300 km/s. This is due 
to the fact that the differences between 77 and 77 , which 
are central to the target dependence of the rate, rapidly 
vanish at Wmin 300 km/s (see Fig. 1). The target- 
independent nature of this region is not specific to mag¬ 
netic DM and occurs whenever the SHM is assumed, at 
least with l/u^ and ^"-dependent terms in the differen¬ 
tial cross section and n ^ 0. This is because all velocity 
Integrals arising from terms going as u" with n ^ 0 in 
the differential cross section have very similar phases at 
ah Vmin values, i.e. they are all comparable to fj In Fig. 1. 
The target-dependent effects addressed in this paper thus 
rely on having both a 1 /w^ term and a u” term, n ^ 0 , 
in the differential cross section. 

At sufficiently small values of Umin the rate Is always 
dominated by r (i.e. / ~ 1 and / ~ 0 ), as shown in 
Fig. 2. This is due to the 1/Umm factor appearing in 
Eq. ( 6 ). Therefore in the small Umin limit one can disre¬ 
gard the contribution of r and correctly assume tmax and 
tmin coincide with the Tmax and Tmin of fj shown In Fig. 1. 
This explains why tmax in Fig. 3 occurs In May at small 
1 /min values regardless of the target. 

Assuming at least one target Isotope has a non-zero 
nuclear magnetic moment, the dipole-dipole part of the 
Interaction becomes dominant, and thus r > f, at large 
values of Umin- This is due to the fact that the spin- 
independent charge form factor decreases faster than the 
magnetic form factor. Fig. 2 confirms that for the ele¬ 
ments and DM masses considered, there is a Umin value 
above which r dominates and below which f dominates. 
In Fig. 3 this corresponds to the time variation of the 
rate being determined by 77 or fj, respectively. For ger¬ 
manium, this switch occurs at large Umin values because 
of its small average magnetic moment. How and where 
this switch in Umin occurs determine the main features of 
^max and tmin ffi Fig- 3. 

For each element, the features in Fig. 2 move to smaller 
Vmin values as the DM particle mass increases. This is in 
part because the Umin value corresponding to a particu¬ 
lar Ar decreases, but also because the 1/fJ-T 
factors in Eq. ( 6 ) decrease. Notice that, as to increases, 
the Umin value above which r becomes the dominant term 
in the rate may fall below 300 km/s, leading to the ap¬ 
pearance of a feature In Fig. 3. This happens with xenon 
when TO goes from 100 GeV to 1 TeV. 

We emphasize that the interplay between 77 and fj does 
not only affect observables associated with the modula¬ 
tion of the rate, such as tmax and tmin, but also the extent 
to which the standard approximation of the modulation 
given In Eq. (1) holds. Fig. 4 shows that the difference be¬ 
tween tmin and tmin = ^max ~ 6 months is target and DM 
particle mass dependent, and can be large, e.g. tmin —tmin 
for TO = 100 GeV could be as large as ±45 days. This 
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FIG. 3. Time of maximum fmax {left) and minimum fmin {right) of the scattering rate as functions of Vmin for a 100 GeV (solid 
lines) and 1 TeV (dashed lines) magnetic dipole DM particle scattering elastically off different targets, assuming the SHM. For 
germanium, the two lines would overlap (only dashed is shown). Also shown are the low energy thresholds for several current 
direct detection experiments, for m = 100 GeV (for larger m the thresholds shift to lower Vmin values). 
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FIG. 4. tmin — tmin, with tmin the time six month apart from 
tmax. See Fig. 3 for details. 


implies that higher order terms in the Fourier expansion 
of the rate beyond Eq. (1) cannot be neglected. 

To illustrate how important the target element depen¬ 
dence of the rate modulation can be, consider the sig¬ 
nal due to a 100 GeV DM particle being detected with 
both xenon and fluorine near the present LUX and PICO 
thresholds. Were the modulation due solely to y or fj, the 
two experiments should observe nearly the same value of 
fmax, see Fig. 1. Instead, due to the target-dependent 
interplay of rj and fj, the tmax observed with the two tar¬ 
get elements could differ by more than four months and 
the modulation in xenon would be better described by 
Eq. (1) than the modulation in fluorine. 


As we already mentioned, in order to observe the 
target-dependent effects described so far, it is essential 
that the experimental threshold in nmin; which depends 
on the threshold in E^, the DM particle mass and the 
scattering kinematics, is below 300 km/s. Fig. 3 shows 
that m = 100 GeV is already large enough with present 
thresholds to observe this effect. For lower m the ef¬ 
fect will only be present with the light targets, for elastic 
scattering. 

Should DM scatter inelastically off nuclei, the scat¬ 
tering kinematics would be different from that of elas¬ 
tic scattering. Inelastic scattering [51, 52] can happen 
if there are at least two almost degenerate DM par¬ 
ticles with masses m and m + S {S ^ m). If the 
particle with mass m scatters into the m + S particle, 
i^min = {{tutEyi/ fir) + S\/^/2mTE^. In particular, if 
(5 < 0 (exothermic scattering [52]), the nmin value cor¬ 
responding to given Ur and m can be much smaller than 
in the case of elastic scattering. 

All the effects we have described here rely on having 
a DM-nucleus differential cross section with a particular 
V dependence. The issue remains of how such a cross 
section could be identified experimentally. We believe 
that this would require observing an annual modulation 
in at least two experiments with different target materi¬ 
als. If the velocity and target dependence in the differ¬ 
ential cross section factorize, the observables associated 
with the modulation as functions of i;niin would be in¬ 
dependent of the target element, for any DM distribu¬ 
tion. However, experiments do not measure their signal 
in Vniin) but in energy, and the values of m and <5 enter¬ 
ing the UR-Wniin relation are not known a priori. This 
problem could be overcome by comparing observables of 
the modulation, like tmax and tmin, of at least two experi- 
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merits employing different target materials, and trying to 
find values of m and 6 that reconcile the differences be¬ 
tween observed modulations as functions of Umin- Should 
there exist no i?R-Ujnin relation that would make the 
modulations as functions of Umin compatible across exper¬ 
iments, one may infer the differential cross section con¬ 
tains a non-factorizable velocity and target dependence. 


CONCLUSIONS 

It is usually assumed that the modulation of the ex¬ 
pected differential rate in direct DM detection experi¬ 
ments, expressed as a function of Umin (the minimum 
DM speed necessary to impart a certain recoil energy to 
a target nucleus), does not depend on the target. We 
have shown instead that experiments employing differ¬ 
ent target materials could observe an entirely different 
annual modulation of their differential rate as a function 
of i^min- This would be a signature of DM interactions 
with more than one velocity-dependent term in the scat¬ 
tering cross section, in particular terms proportional to 
1/u^ and v'^ with n ^ 0. In order to identify experimen¬ 
tally this type of cross section, we believe at least two 
experiments employing different target materials should 
observe an annual modulation. Should no UlR-Umin re¬ 
lation be found that reconciles the modulated signals as 
functions of Umin, one may infer the differential cross sec¬ 
tion contains a non-factorizable velocity and target de¬ 
pendence regardless of the DM distribution. 

As an example, we have shown explicitly the target 
dependence of the time of maximum t„iax and minimum 
fmin of the rate for a 100 GeV and a 1 TeV magnetic 
dipole DM scattering elastically assuming the SHM. We 
found that the values of tmax observed with e.g. xenon 
and fluorine close to the present LUX and PICO thresh¬ 
olds could disagree by as much as four months (see 
Fig. 3), and the modulation in xenon could be better 
described by the sinusoidal time dependence usually as¬ 
sumed than that in fluorine. 


ACKNOWLEDGMENTS 

E.D.N. and G.G. acknowledge partial support from 
the Department of Energy under Award Number DE- 
SC0009937. 


[1] A. Drukier, K. Freese, and D. Spergel, Phys.Rev. D33, 
3495 (1986b 

[2] S. K. Lee, M. Lisanti, and B. R. Safdi, JCAP 1311, 033 
(2013), arXiv: 1307.5323 [hep-ph]. 

[3] P. Smith and J. Lewin, Pliys.Rept. 187, 203 (1990). 


[4] M. S. Alenazi and P. Gondolo, Phys.Rev. D74, 083518 
(2006), arXiv:astro-ph/0608390 [astro-ph]. 

[5] S. K. Lee, M. Lisanti, A. H. G. Peter, and B. R. Safdi, 
Phys.Rev.Lett. 112, 011301 (2014), arXiv: 1308.1953 
[astro-ph. GO]. 

[6] M. Pospelov and T. ter Veldhnis, Phys.Lett. B480, 181 
(2000), arXiv:hep-ph/0003010 [hep-ph]. 

[7] K. Sigurdson, M. Doran, A. Knrylov, R. R. Caldwell, 
and M. Kamionkowski, Phys.Rev. D70, 083501 (2004), 
arXiv:astro-ph / 0406355 [astro-ph] . 

[8] V. Barger, W.-Y. Kenng, and D. Marfatia, Phys.Lett. 
B696, 74 (2011), arXiv: 1007.4345 [hep-ph]. 

[9] S. Chang, N. Weiner, and I. Yavin, Phys.Rev. D82, 
125011 (2010), arXiv: 1007.4200 [hep-ph]. 

[10] W. S. Cho, J.-H. Hnh, I.-W. Kim, J. E. Kim, and 
B. Kyae, Phys.Lett. B687, 6 (2010), arXiv: 1001.0579 
[hep-ph]. 

[11] J. H. Heo, Phys.Lett. B693, 255 (2010), arXiv:0901.3815 
[hep-ph]. 

[12] S. Gardner, Phys.Rev. D79, 055007 (2009), 

arXiv:0811.0967 [hep-ph]. 

[13] E. Masso, S. Mohanty, and S. Rao, Phys.Rev. D80, 
036009 (2009), arXiv:0906.1979 [hep-ph]. 

[14] T. Banks, J.-F. Fortin, and S. Thomas, (2010), 
arXiv: 1007.5515 [hep-ph]. 

[15] J.-F. Fortin and T. M. Tait, Phys.Rev. D85, 063506 
(2012), arXiv:1103.3289 [hep-ph]. 

[16] H. An, S.-L. Chen, R. N. Mohapatra, S. Nussi- 
nov, and Y. Zhang, Phys.Rev. D82, 023533 (2010), 
arXiv: 1004.3296 [hep-ph]. 

[17] K. Kumar, A. Menon, and T. M. Tait, JHEP 1202, 131 
(2012), arXiv:1111.2336 [hep-ph]. 

[18] V. Barger, W.-Y. Kenng, D. Marfatia, and P.-Y. Tseng, 
Phys.Lett. B717, 219 (2012), arXiv: 1206.0640 [hep-ph]. 

[19] E. Del Nobile, C. Kouvaris, P. Panel, F. Sannino, and 
J. Virkajarvi, JCAP 1208, 010 (2012), arXiv: 1203.6652 
[hep-ph]. 

[20] J. M. Cline, Z. Lin, and W. Xue, Phys.Rev. D85, 101302 
(2012), arXiv:1201.4858 [hep-ph]. 

[21] N. Weiner and 1. Yavin, Phys.Rev. D86, 075021 (2012), 
arXiv: 1206.2910 [hep-ph]. 

[22] S. Tulin, H.-B. Yu, and K. M. Zurek, Phys.Rev. D87, 
036011 (2013), arXiv:1208.0009 [hep-ph]. 

[23] J. M. Cline, A. R. Frey, and G. D. Moore, Phys.Rev. 
D86, 115013 (2012), arXiv: 1208.2685 [hep-ph]. 

[24] J. H. Heo and C. Kim, Phys.Rev. D87, 013007 (2013), 
arXiv: 1207.1341 [astro-ph.HE]. 

[25] E. Del Nobile, G. Gelmini, P. Gondolo, and J.-H. Huh, 
JCAP 1310, 048 (2013), arXiv:1306.5273 [hep-ph]. 

[26] G. Barello, S. Chang, and C. A. Newby, Phys.Rev. D90, 
094027 (2014), arXiv:1409.0536 [hep-ph]. 

[27] I. Lopes, K. Kadota, and J. Silk, Astrophys.J.Lett. 780, 
L15 (2014), arXiv:1310.0673 [astro-ph.SR]. 

[28] W.-Y. Keung, 1. Low, and G. Shaughnessy, Phys.Rev. 
D82, 115019 (2010), arXiv:1010.1774 [hep-ph]. 

[29] M. I. Gresham and K. M. Zurek, Phys.Rev. D89, 016017 
(2014), arXiv:1311.2082 [hep-ph]. 

[30] E. Del Nobile, G. B. Gelmini, P. Gondolo, and J.-H. 
Huh, JCAP 1406, 002 (2014), arXiv: 1401.4508 [hep-ph]. 

[31] M. 1. Gresham and K. M. Zurek, (2014), arXiv:1401.3739 
[hep-ph]. 

[32] C. M. Ho and R. J. Scherrer, Phys.Lett. B722, 341 
(2013), arXiv:1211.0503 [hep-ph]. 



6 


[33] A. L. Fitzpatrick and K. M. Zurek, Phys.Rev. D82, 
075004 (2010), arXiv: 1007.5325 [hep-ph], 

[34] M. T. Frandsen, F. Kahlhoefer, C. McCabe, S. Sarkar, 
and K. Schmidt-Hoberg, JCAP 1307, 023 (2013), 
arXiv:1304.6066 [hep-ph], 

[35] Y. Gao, C. M. Ho, and R. J. Scherrer, Phys.Rev. D89, 
045006 (2014), arXiv:1311.5630 [hep-ph]. 

[36] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, 
T. M. Tail, et al, Phys.Rev. D82, 116010 (2010), 
arXiv:1008.1783 [hep-ph]. 

[37] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, 
T. M. Tail, et al, Phys.Lett. B695, 185 (2011), 
arXiv:1005.1286 ]hep-ph]. 

[38] J.-M. Zheng, Z.-H. Yu, J.-W. Shao, X.-J. Bi, Z. Li, et al, 
Nucl.Phys. B854, 350 (2012), arXiv:1012.2022 [hep-ph]. 

[39] Z.-L. Liang and Y.-L. Wu, Phys.Rev. D89, 013010 
(2014), arXiv:1308.5897 [hep-ph]. 

[40] R. Gatena and P. Gondolo, (2014), arXiv:1405.2637 [hep- 
ph]. 

[41] R. Catena, JCAP 1407, 055 (2014), arXiv:1406.0524 
[hep-ph]. 

[42] A. L. Fitzpatrick, W. Haxton, E. Katz, N. Lubbers, and 
Y. Xu, JCAP 1302, 004 (2013), arXiv:1203.3542 [hep- 
ph]. 

[43] A. L. Fitzpatrick, W. Haxton, E. Katz, N. Lubbers, and 
Y. Xu, (2012), arXiv:1211.2818 [hep-ph]. 

[44] J. Kumar and D. Marfatia, Phys.Rev. D88, 014035 
(2013), arXiv:1305.1611 [hep-ph]. 

[45] M. Cirelli, E. Del Nobile, and P. Panel, JCAP 1310, 019 
(2013), arXiv:1307.5955 [hep-ph]. 

[46] G. B. Gelmini, A. Georgescu, and J.-H. Huh, JGAP 
1407, 028 (2014), arXiv:1404.7484 [hep-ph]. 

[47] D. Akerib et al. (LUX Gollaboration), Phys.Rev.Lett. 
112, 091303 (2014), arXiv:1310.8214 [astro-ph.CO]. 

[48] R. Agnese et al. (SuperCDMS Collaboration), (2014), 
arXiv:1402.7137 [hep-ex]. 

[49] R. Bernabei et al, (2013), arXiv: 1308.5109 [astro- 
ph.GA]. 

[50] C. Amole et al. (PICO), (2015), arXiv: 1503.00008 [astro- 
ph.CO]. 

[51] D. Tucker-Smith and N. Weiner, Phys.Rev. D64, 043502 
(2001), arXiv:hep-ph/0101138 [hep-ph]. 

[52] P. W. Graham, R. Harnik, S. Rajendran, and 

P. Saraswat, Phys.Rev. D82, 063512 (2010), 

arXiv:1004.0937 [hep-ph]. 


